Harris sheet solution for magnetized quantum 

plasmas 

F. Haas* 

Universidade do Vale do Rio dos Sinos - UNISINOS 
Unidade de Exatas e Tecnologicas 
Av. Unisinos, 950 
93022-000 Sao Leopoldo, RS, Brazil 

February 2, 2008 



Abstract 

We construct an infinite family of one-dimensional equilibrium so- 
lutions for purely magnetized quantum plasmas described by the quan- 
tum hydrodynamic model. The equilibria depends on the solution of 
a third-order ordinary differential equation, which is written in terms 
of two free functions. One of these free functions is associated to the 
magnetic field configuration, while the other is specified by an equa- 
tion of state. The case of a Harris sheet type magnetic field, together 
with an isothermal distribution, is treated in detail. In contrast to the 
classical Harris sheet solution, the quantum case exhibits an oscilla- 
tory pattern for the density. 



1 Introduction 

Quantum plasmas have attracted renewed attention in the last years, due e.g. 
to the relevance of quantum effects in ultra-small semiconductor devices PQ, 
dense plasmas [2] and very intense laser plasmas |3J. The most recent devel- 
opments in collective effects in quantum plasmas comprises wave propagation 
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in dusty plasmas [I]-[H], soliton and vortex solutions PTQ|, [TTj . shielding effects 
[T2"l H3], modulational instabilities [H] and spin effects [15]. Most of these 
works have been made using the hydrodynamic model for quantum plasmas 
[16]- [20], in contrast to more traditional approaches based on kinetic descrip- 
tions [21]. Microscopic descriptions like coupled Schrodinger equations or 
Wigner function approaches are more expensive, both numerically and ana- 
lytically, specially if magnetic fields are allowed. For a general review on the 
available quantum plasma models, see [22]. 

The electrostatic fluid model for quantum plasmas have been recently 
extended to incorporate magnetic fields |23j. The new quantum hydrody- 
namic model was derived taking the first two moments of the electromag- 
netic Wigner equation, which is the quantum counterpart of the correspond- 
ing Vlasov equation, and assuming a closure condition p = p(n). In other 
words, the procedure is formally the same as for classical plasma fluid de- 
scriptions, while now the starting point is the Wigner-Maxwell and not the 
Vlasov-Maxwell system. The electromagnetic quantum fluid model has been 
already used for the analysis of shear Alfven modes in ultra-cold quantum 
magnetoplasmas [21], the description of drift modes in nonuniform quantum 
magnetoplasmas [25] and of shear electromagnetic waves in electron-positron 
plasmas [26] . Instead of the discussion of wave propagation in quantum plas- 
mas, the aim of this letter is the analysis of some simple quantum magneto- 
static equilibria resembling the well known Harris profile for classical plasma 



2 Quantum magnetoplasma equilibria 

For a one-component quantum plasma, the electromagnetic quantum fluid 
equations reads [23] 



All the symbols in eqs. (JT]]2]) have their conventional meaning and the system 
is supplemented by Maxwell equations. Only electrons are considered, the 
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ions being described by a convenient immobile background. Notice the extra 
dispersive term, proportional to h 2 , at the moment transport equation. This 
Bohm potential term has profound consequences on the structure of the 
equilibrium solutions, as we shall see in the following. 

Specifically, consider a purely magnetic one-dimensional class of time- 
independent solutions characterized by zero electric field and 

B = B y (x)y + B z [x)z , 

n = n(x) , (3) 
u = u y (x)y + u z {x)z , 
p = p(n) . 

The magnetic field can be given in terms of a vector potential A = A y (x)y + 
A z (x)z, so that By = —dA z /dx and B z = dA y /dx. Notice that the fluid 
model is suitable for the search for static quantum equilibria since the kinetic 
(Wigner) equation is not satisfied by arbitrary functions of the invariants of 
motion as for Vlasov plasmas. Therefore we are not allowed to use Jeans 
theorem for the construction of equilibria. 

Neutrality is enforced by an appropriate immobile ionic background de- 
scribed by an ionic density rii{x). Therefore, Poisson's equation can be ig- 
nored. Now inserting the form ([3]) into Ampere's law and the quantum fluid 
equations gives 

cPA,. 



dx 2 
d 2 A z 

dx 2 



efi nu y , (4) 
e/j, nu z , (5) 



dp . dAy dA z 

— = - en {u^- + u z ——) 

dx ' dx dx 

h 2 n d ( d 2 y^n/dx 2 \ 
2m dx \ y/n J 

As in the classical situation [28J, it is useful to restrict to the cases 
where the magnetic field is indirectly defined through a pseudo-potential 
V = V(A y , A z ) for which 

1 dV 



nu y 



(7) 



e/xo dAy 
1 dV 

nu z = — — , (8) 

efjbo dA z 
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so that (jMS]) is transformed into a two-dimensional autonomous Hamiltonian 

system, 



d 2 A 



dV 

dAy 

dV 



(9) 



dx 2 
d 2 A z 



(10) 



dx 2 



In this system, x plays the role of time, while the components of the vector 
potential play the role of spatial coordinates. After specifying the pseudo- 
potential V and solving eqs. (l9l fT0|) . we regain the magnetic field using B = 
V x A. The current then follows from eqs. (J7HEJ)- 

The choice expressed at eqs. (JIIISJ) imposes a restriction on the classes 
of equilibria, since not all density and velocity fields can be cast in this 
potential form. However, introducing the pseudo-potential V has at least two 
advantages. First, we can learn from Hamiltonian dynamics how to design 
specific pseudo-potentials V in order to obtain special classes of magnetic 
fields. For instance, periodic magnetic fields can be easily obtained from well 
known potentials associated to periodic solutions. Second, the formalism 
becomes more compact in terms of the function V. 

In terms of V, the balance eq. ([6]) reads 



It can be shown that, apart from an irrelevant numerical constant, the 
pseudo-potential V is directly related to magnetic pressure, V = — |B| 2 /2, 
showing that the left-hand side of eq. ( TTTT) refers to the usual (classical) 
pressure balance equation. The right-hand side, however, has a pure quan- 
tum nature. Not only there must be a balance between kinetic and magnetic 
pressures, since the quantum pressure arising from the Bohm potential term 
has to be taken into account. This quantum pressure manifests e.g. in 
the dispersion of wave-packets in standard quantum mechanics. In plasmas, 
the quantum pressure is responsible for subtle effects like in the case of the 
quantum two-stream instability, where the instability is magnified for small 
wave-numbers and suppressed for large wave- numbers [TH [20] . 

In the quantum case where h ^ 0, eq. ( TTTT) is a third-order ODE for the 
density. It is useful to express this equation in terms of a variable a = \fn. 




(11) 
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Taking into account the equation of state p = p(n) and denning a new 
function V(x) = V(A y (x), A z (x)), we get 

aa'" - a'a" + f(a)a' + g(x) = 0, (12) 

where the prime denotes differentiation with respect to x and we have intro- 
duced the quantities 

Amadp 2 
f(a) = -^-^(„ = a), (13) 

9{x) = ^rL- (14) 

The strategy to derive the solutions is now clear. Choosing a pseudo-potential 
V(A y , A z ) and then solving the Hamiltonian system flUl fTUj) for the vector 
potential, we determine simultaneously the magnetic field and V. Quantum 
effects manifests in the equation for the density, eq. (TT2l . which also deserves 
the function of state p = p(n). 

Another legitimate interpretation of the balance equation (fl2i) is to first 
specify the particle density n and the magnetic pressure — V/no and then 
solving for the kinetic pressure. This would give an equation of state with a 
quantum correction. However, in most applications, one supposes a certain 
equation of state and then proceeds to the calculation of the density and 
velocity fields. This will be our preferred approach in what follows. In the 
next section, we consider in detail the case of Harris sheet magnetic fields. 



3 Quantum Harris sheet 

Exactly as for the classical Harris solution, suppose a isothermal plasma, 
p = nn B T, and a pseudo-potential function 

V=^eM^j), (15) 

where L is a characteristic length and is a (constant) magnetic field 
reference value. The Hamiltonian system flUl fTU]) is then 

^ = °' lt^ = -— eM BlL h (16) 
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If we further take the boundary conditions A z (x = 0) = [dA z j dx)[x = 0) = 
0, we easily solve (TT5|) to get 



A y = A y0 + B x , A z = —B^L In cosh(x/ L) , (17) 

where A y o and Bq are integration constants. The magnetic field following 
from this vector potential characterizes the well-known Harris sheet solution, 

B y = B^ tanh(x/L) , B Z = B , (18) 

also allowing for a superimposed homogeneous magnetic field. 
In addition, the velocity field follows from 



u y = 0, u z = B °° sech 2 (^-) . (19) 
efj, Q nL L 

Notice that any departure from the classical density solution would imply 
further changes in the velocity field. 

To derive the density we have to solve the third-order ODE eq. (fl2]) . 
constructed in terms of the functions f(a) and g(x) at (I13til4j) . Using the 
isothermal equation of state, the form ( TloT) for the pseudo-potential V and 
the Harris sheet solution, we get 

f[a) = — 2 — a, (20) 

g(x) = -^ sech2 (|) tanh (f)- (21) 

Adopting the dimensionless variables 

a = a/y/no, X = x/L, (22) 
where n is some ambient density such that 

n K B T = 3° , (23) 

eq. ([HD is finally expressed as 

d?a da d 2 a a da 1 i2aa , n 

a dx* ~ dxdx* ~lPdx = lp sech xtanh x ' (24) 
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in terms of a new dimensionless parameter 



where V a = -B 00 /(/io??Wo) 1//2 is the Alfven velocity. 

The parameter H is a measure of the relevance of the quantum effects. 
It is essentially the ratio of the scaled Planck constant h to the action of 
a particle of mass m travelling with the Alfven velocity and confined in a 
length L related to the thickness of the sheet. The larger the ambient density 
n and the smaller the characteristic length L or the characteristic magnetic 
field Boo, the larger are the quantum effects. 

In order to understand the role of the quantum terms, we may investigate 
dMD with 

c\n (fin 

«(X = 0) = 1, ^(X = 0) = 0, |£(X = 0) = -1, (26) 

which reproduces the boundary conditions for the classical Harris sheet, when 
a = sechX. With the choice fl26l) . eq. ( 1211) integrated once gives 




-sech 2 x) . (27) 



In the ultra-quantum limit H —>■ oo, the left-hand side of fT27]) vanishes. In 
this situation and using the prescribed boundary conditions, the solution is 

a = cos X. (28) 

This imply a qualitative change (from localized to oscillatory) on the solution 
due to quantum effects. In order to further investigate this tendency, we show 
the numerical solution for ff27]) with the appropriate boundary conditions for 
a few values of H. This is shown in the figs. 1 and 2, where increasingly 
oscillatory solutions are shown, according to H = 1 or H = 5 respectively. 
On the opposite case, fl2TI) shows that when H — >• we regain the classical 
Harris solution, a = sechX. 

Another interesting possibility is an equation of state for an ultra-cold 
Fermi gas, 

5n ' 
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where Tp is the Fermi temperature and n is the ambient density. Proceeding 
exactly as before and assuming 



n K B T F = ^ , (30) 

we obtain 

d 3 a da d 2 a a 7 ^ 3 da 1 2 

a dx 3 ~ dxdx 2 ~lPdx = lP sech x tanh x ' (31) 

where a and X are defined in ff22l) and H in fl25l) . Similar oscillatory behaviour 
is also found for nonzero H and suitable boundary conditions. 



4 Summary 

Equation (fl2|) describes a whole class of quantum equilibria for magnetoplas- 
mas. The particle density compatible with a tanh magnetic field shows an 
increasingly oscillatory pattern, in comparison to the classical system associ- 
ated to a localized sech 2 solution. Other classes of equilibria can be built for 
different choices of pseudo-potentials V = V(A y , A z ) and equations of state 
p = p(n). The ideas in the present formulation may be a starting point for 
magnetic equilibria relevant for dense astrophysical objects like white dwarfs. 
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